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Abstract 

Logistic growth on a static heterogenous substrate is studied both above and below the drift- 
induced derealization transition. Using stochastic, agent-based simulations the derealization of 
the highest eigenfunction is connected with the large N limit of the stochastic theory, as the local- 
ization length of the deterministic theory controls the divergence of the spatial correlation length 
at the transition. Any finite colony made of discrete agents is washed away from a heterogeneity 
with compact support in the presence of strong wind, thus the transition belongs to the directed 
percolation universality class. Some of the difficulties in the analysis of the extinction transition 
in the presence of a localized active state are discussed. 
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The effect of drift on inhomogeneous systems that exhibit growth and propagation has 
attracted much interest in the last decade Q, [3, [3, Q, [f], 3, H Isl, [ill 11 1- When the time 
evolution of a system is governed by a real symmetric evolution operator it may support 
both extended and localized eigenstates. The eigenstates of a quantum particle in a sin- 
gle potential well, for example, are either localized inside the wall or extended above some 
threshold energy. In the presence of drift, or other non-Hermitian perturbation [ij, the sys- 
tem undergoes a phase transition where localized wavefunctions become extended, and the 
corresponding eigenvalues migrate from the real axis to the complex plane. This transition 
was first analyzed by Hatano and Nelson [2] in the context of flux lines in high T c super- 
conductors with columnar defects subjected to a tilted external magnetic field. Since then, 
many authors have considered this transition in different fields, e.g. hydrodynamics jjj, 



random lasers 



4|, and quantum dots 5] among many others. 



Of particular interest, both theoretically 



3, El and experimentally 
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ll|] , is the 



derealization transition for bacterial colonies on a heterogeneous substrate in the presence 
of drift. A logistic growth of a motile population on a Id static spatially heterogenous 
substrate is described by: 



dc(x, t) 



DV 2 c(x,t) +v 



dc(x, t) 



+ a(x)c(x, t) — c 2 (x, t). 



(1) 



dt dx 

In the absence of drift term (v = 0) and for a homogenous environment (a = ao = cr — fi, 
where ao, the difference between the birth rate a and the death rate /1, is independent 



of spatial location 



121 ]) one gets the celebrated Fisher-Kolomogorov-Petrovsky-Piscounov 



equation (FKPP), a generic description of an invasion of a stable state (c* = ao) into an 
unstable one c* = 0. In the asymptotic long-time limit this system supports a front that 



travels with constant speed vf = 1\jDa§. In the homogenous case the eigenstates of the 
linearized evolution operator 



C = DV 2 c(x, t) + v dC ^ ^ + a c(x, t). 

ox 



(2) 



are extended sinusoidal functions and the drift corresponds to a simple Galilean transfor- 
mation. 

Things change when translational invariance is broken, i.e., in the presence of spatial 
inhomogeneity. Two main types of heterogenous growth are considered in the literature 
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ll|: A "single oasis" case, where the growth rate is larger on a spatial domain, and 



the disordered case, where a(x) = ao + Sa(x), 5a being taken from some random distribution 
with zero mean. In both cases the spectrum of the linear evolution operator admits localized 
wavef unctions; if {0°(x),r„} is the set of eigenfunctions and the corresponding eigenvalues 
of £(v = 0), at least some eigenstates in the tail of the spectrum (or all the states for a 
disordered potential below 2d) are exponentially localized. The effect of small drift on a 
localized eigenstates is trivial: 

= e ra/2D 0°(x) , K = T° n -^. (3) 

This "gauge invariance" breaks down at v c n = 2D/£ n , where £ n is the localization length 
of the nth eigenstate. Above v c n the eigenstate delocalizes and the boundary conditions 
begin to play an important role: e.g., for periodic boundary conditions the eigenvalues that 
correspond to delocalized eigenstates become complex 13j. The spectrum then takes the 
form of a "bubble" in the complex plane, where the localized eigenstates correspond to 
the spectral points in the tail, since the localization length in the center of the band is 
larger. A non-Hermitian "mobility edge" appears between the two regimes. Increasing v 

points, and at the 



15j ensures that the 



even more, the bubble spreads and captures more and more spectra 
end the ground state also delocalizes. The Perron- Frobenius Theorem 
highest eigenstate stays on the real line, and the derealization transition is identified by the 

nn 

breakdown of the trivial gauge (Eq. [3]) and the vanishing of the spectral gap |6|, 1 131 ] - 

Figure 1 shows some examples of the spectrum of C, together with a sketch of the phase 
diagram, for the single oasis scenario. In the absence of drift there is a single localized state 
at the right edge of the spectrum (if the oasis is large, a few localized states exist), followed 
by a continuum of states that correspond to extended eigenfunctions. Even a small drift 
is enough to push the delocalized eigenvalues to the complex plane, but the localized state 
only develops a slight asymmetry with almost no effect on the eigenvalue. Only for high 
enough drift does the highest eigenstate delocalize and the gap disappear. A change of a 
corresponds to a rigid shift of the whole spectrum along the real line. Thus, three regimes 
exist in the drift-proliferation parameter space: the extinction region, where the real part 
of all the eigenvalues is negative; the localized region, where only the localized states admit 
positive growth rate; and the proliferation regime, where both localized and extended states 
may grow. Above Vq only the extinction and the proliferation regions exist. 

The above discussion is, however, too naive. Bacterial systems are not deterministic, and 
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FIG. 1: A sketch of the "phase diagram" for a single oasis in the death-rate - drift space. Region 
(1) is the extinction region where the highest eigenvalue is negative. In region (2) only the localized 
state admits an eigenvalues with positive real part, and in regions (3a) and (3b) extended eigenstates 
become "active". In the right side of the figure the highest eigenstate and the spectrum in the 
complex plane are plotted for the cases of no drift (v = 0), small drift (v = 0.1) and large drift 
(v = 0.5) for a single oasis. 



are composed of discrete objects that may die, reproduce or migrate with some probability 
that depends on the local environmental conditions. The bacterial population at a certain 
point is not a deterministically varying continuous quantity, like c(x), but a discrete number 
that undergoes stochastic processes, e.g., A — > 2 A, A + A — > 0, A ^ etc. Like in many 
other branches of science, the deterministic dynamics is an approximate description of the 
system that becomes exact where the effect of stochasticity vanishes. In the case considered 
here the demographic stochasticity becomes negligible when the density of agents is large, 
since the relative fluctuations scale with 1/yN. Technically, the exact stochastic Master 
equation is replaced by a deterministic description using the Kramers- Moyal expansion, 
or more rigorously by van-Kempen's Q expansion and related methods 



Lebowitz 
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Joo and 



16( | have already pointed out that in the limit of large N one should expect a 
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population density distribution that follows the spatial features of the active eigenstates, 
i.e., the eigenstates for which Re(T) > 0. Here, on the other hand, we want to discuss 
the effects of spatial heterogeneity and drift for a dilute system; i.e., close to the extinction 
transition. In that case the 1 / y/N expansions are invalid and so we resort to numerical 
simulations. 

Let us first present some general considerations. Grassberger and Janssen [l^] suggested 
long ago that the extinction transition to a single absorbing state on a homogenous sub- 
strate belongs (in the absence of special additional symmetries) to the directed percolation 
(DP) equivalence class, independent of the microscopic details of the stochastic process. DP 
is a continuous transition and the correlation length and correlation times diverge at the 
transition point with their characteristic exponents (see for a general review). On a 
homogenous substrate the correlation length is the only length scale of the problem. On a 
static heterogenous substrate, on the other hand, another length scale appears - the local- 
ization length. How do these two quantities relate to each other? What are the properties of 
the stochastic extinction transition below and above the deterministic derealization transi- 
tion? In what sense is Eq. (CQ) a deterministic limit of a stochastic process when the effect 
of stochasticity is important; i.e., close to the extinction transition? 

Recently, this last question has been addressed for the transition on a homogenous sub- 



strate 



211 ] . It turns out that the transition is always in the DP equivalence class, but the 



carrying capacity of the system, N, determines the location of the transition and, more im- 
portant, the width of the transition zone. In the deterministic theory the correlation length 
is zero both below and above the transition (any initial density fluctuation simply decays 
exponentially to the stable state and its spread during this process is negligible). The spatial 
correlation length for the stochastic process satisfies £j_ ~ A~ U± /N K , where A is the distance 
from the transition. If the system parameters are such that its deterministic analogue is at 
the transition point, then A (AT) ~ N~ T and A vanishes at the deterministic limit. Under 
these conditions £j_ ~ iV 7 , where 7 = tv± — k. This implies that for any finite A for large 
enough N the correlation length shrinks to zero and the deterministic description holds, 
whereas for any finite N, for small enough A the system enters the transition zone and the 
deterministic description collapses. Both k and r depend on the deterministic features of 



the model as explained in 
at the transition point. 



21|, but in any case 7 > so the deterministic limit never exists 
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In order to simulate the heterogenous system in the large N limit, an individual based 
model that allows for an accurate determination of the transition point in the N — >• oo limit 
is used. We consider a logistic growth process on a one dimensional lattice with periodic 
boundary conditions; Euler integration is used with small, but finite, At. The number if 
agents at the 2-th lattice site is an integer rij, and each cycle of the Monte-Carlo simulation 
involves two consecutive steps. The first step is the reaction: each of the agents at the site 
produces an offspring with probability (ctq + &7j)(l — ni/No)At, and dies with probability 
\i At. In the second, diffusion step, any agent is selected for migration with probability 
2\At, then chooses its destination - to the left with probability qi, = (1 + v)/2 or to the 
right with probability 1 — q^. To avoid artificial drift as a result of the sequential update 
of lattice sites, parallel update was used; is updated only after the diffusion cycle is 
completed. 

In the linearized deterministic limit this model corresponds to an L dimensional map, 
where L is the number of sites. This map is given by the multiplication of the reaction 
matrix, Rij = Sij[l + At(oo + 8(Ti — /i)], by the diffusion matrix that takes the form (up 
to the boundary conditions) D t j = 5^(1 — 2% At) + + u)At + <5i+i,jx(l — v)At. 

Diagonalizing the product DR one finds the highest eigenvalue T and the corresponding 
eigenvector, <j) l ; adding another death process, where each particle in the MC simulation is 
selected to die after any cycle with probability 1/ 

tildeGammao, ensures that the system is exactly at the transition point for iVo — > oo. 
In different words, the agent-based system is simulated with a parameter set that ensures 
To = in the deterministic limit. 

Clearly, a system with finite carrying capacity N is always closer to extinction than the 
deterministic system when all other parameters are equal. This implies that, scaling the 
parameters as described above and increasing N , the system is always in the extinction 
phase and reaches the transition exactly at Nq = oo. In Fig. EE the correlation length, 
is plotted vs. iVo on a log-log scale and reveals the real meaning of the deterministic der- 
ealization transition: below v c , i.e., when 0o is localized, the correlation length associated 
with the stochastic process first grows and then saturates to the deterministic localization 
length £/[l — v£/(2D)]. This demonstrates the fact that the state that becomes active at 
the transition is localized and the correlation length of the stochastic process can not grow 
beyond this deterministic length. 
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FIG. 2: The log of the transverse correlation length £j_ (in the direction of the drift) vs. log(No) 
for a system of length L = 2000. Here £ = 0.2, do = /i = 1 and the heterogeneity 5ai = 0.26^^/2- 

On the other hand, above v c the correlation length grows unboundedly with iVo- The 
data is consistent with £j_ ~ A*^' 5 , which is (up to logarithmic corrections), the lifetime of 



a well-mixed system at the deterministic transition point [20j. This reflects the fact that 
the delocalized system is not really one dimensional but rather "0 + 1" dimensional, with 
the spatial direction playing the role of time. Life in that system is a result of a drift from 
a source, not of uniform growth, and the lifetime of the colony at the oasis determines the 
spatial extent reached by its decedents. 

For finite iVo and with constant drift velocity v, the system undergoes undergone an 
extinction transition as ao decreases. This may happen either via the localized phase (e.g., 
along Path 1 shown by the arrow in Figure 1), or directly to the delocalized phase (Path 2 
in figure 1). While for A" — > oo the transition happens when i?e(r ) touches zero, for any 
finite Ao the transition takes place when a finite region of the upper part of the spectrum 
is above zero (inset of Fig. 3). For a single oasis (or otherwise when the number of oases is 
finite) all the localized states decay in the long run as a result of demographic stochasticity; 
only when the extended eigenstates are "excited" (their eigenvalues cross to the positive real 
part of the spectrum) will the system be in its active phase. As a result, the scenarios 1 and 
2 can be seen to differ significantly. 

Let us first consider path 2. Intuitively, above v c the colony is carried off the oasis by 
the wind, thus the large-scale properties of the system are identical with a homogenous 
substrate with drift in the thermodynamic limit. More precisely, for finite N the transition 
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occurs when a finite part of the spectrum, made of delocalized states, is already "excited" 
(i.e., Re(T) > for the these eigenstates). Thus, there are two regimes. Deep in the 
extinction phase all states decay, Re(T) < 0. The bacterial density in this regime satisfies 
the deterministic solution c(x,t) = exp(— |ao|£ — (x — Xo — vt) 2 / ADt) / V AitDt , where Xq is 
the nucleation point. The overall occupation of a point, Cr(x) = / c(x,t) dt, is thus a 
monotonically decreasing function of x — xq, with an exponential decay of the tail Ct(x) ~ 



exp(— (x — x )/O) where the localization length £ scales like D/(^Jv 2 + 2_D|a | — v). 

Close to the transition point for finite Nq, on the other hand, many linear states are 
already excited and the growth of the colony is unaffected by the nonlinear competition at 
short times. Only after the characteristic time £y does nonlinearity suppress the growth, 
leading to extinction. Within this growth period the system behaves deterministically and 
a "Fisher front" starts to invade the empty region. As the wind velocity is larger than the 
Fisher velocity above v c {(J, the maximum of Ct(x) moves in the direction of the wind, 
as demonstrated in Figure 3. This second regime vanishes at the deterministic limit; ac- 
cordingly, the detachment of the peak from the nucleation point disappears upon increasing 
N . 

The situation is completely different along path 1. The highest state is now localized, and 
its nonlinear interaction differs substantially from the interaction between extended states. 
If the localization length is finite the oasis region decouples from the rest of the system in the 
thermodynamic limit and the DP dynamics happens in parallel with the zero dimensional 
stochastic process on the oasis. This decoupling, however, is impossible at the bulk DP 
transition, when £j_ diverges 22J. A related issue is the transition in the presence of a 
finite density of randomly distributed oases: below v c a nonuniversal Griffiths phase appears 



between the active and the inactive parameter regions [23j. In the deterministic limit only 
the highest localized state becomes active at the transition, thus the Griffith phase admits 
no deterministic limit, and its width shrinks to zero. These last two observations suggest 
that the deterministic description of the system by means of excited localized eigenstates 
is insufficient, as the convergence of a finite N system to the deterministic limit is a very 
subtle issue, to be addressed in subsequent publication. 
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FIG. 3: Ct{x) vs. x close to the transition {(3 = 0.19, v = 0.7, other parameters identical with 
those used in Fig. 2) for Nq = 3. The effect of the localized heterogeneity at x = is almost unseen. 
As emphasized in the upper inset, the transition takes place when many delocalized modes are also 
excited; as a result the system initially grows and the colony pushed to the right by the wind, 
yielding a hump away from the oasis. This peak disappears for larger (3 and for larger Nq (results 
not shown). The lower inset exemplify the spectrum of the linearized evolution operator when the 
transition takes place along path 1 of Figure 1. 
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